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INVESTIGATION  OF  VELOCITY  FIELD  BEHIND  LOW-FLYING  WING  WITH  F'JSEIAGF 


(Presented  at  Scientific  and  Technical  Conference  of  Lenin  Ship 
Euilding  Institute  (L  K I ) , Cctoter  1969) 


M.  K.  Bochin,  Dept,  of  Hydromechanics;  L.  G.  Tsvetkov,  Dept,  of 
Applied  and  Computer  Mathematics 


In  this  article  we  propose  a method  cf  investigating  the 
velocity  field  behind  a low-flying  wing  with  fuselage.  The  viluc  cf 
the  velocity  field  is  necessary  in  calculating  the  aerodynamic 
characteristics  of  tne  tail  assembly  of  aircraft  which  use  the  screen 
effect. 
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Numerical  Calculation  Method  of  Velocity  Field 

Vie  investigate  the  uniform  movement  of  a wing  of  tinite  span 
with  a fuselage  above  a solid  screen,  lie  sheet  of  solid  vertices 
which  is  formed  twist,  according  to  [1],  into  two  vortex  tip 
filaments,  which  lie  upstream  from  the  wing.  The  intensity  ana 
arrangement  of  the  vortices  are  determined  in  the  main  velocity  fieli 
behind  the  wing.  Therefore,  in  our  study  we  will  ignore  the  stream 
causec  by  the  attached  wing  vortex.  Cccidirate  system  ox'v'z', 
attached  to  the  wing  (Fig.  1),  is  selected  such  that  the  attache! 
vortex  of  the  supporting  wing  lies  within  plane  y'oz';  the  ex*  axis 
lies  downstream  from  the  fuselage;  the  cy ' axis  vertically  abov  it. 
The  fuselage  is  assumed  to  have  the  shape  cf  a circular  finit*- 
cylinder  with  a zero  angle  cf  attack.  Cur  symbols  are:  u0  - velocity 
cf  undisturbed  oncoming  flew  parallel  + c ox*  axis  and  plane  of 
screen;  - intensity  of  vortf  x tip  filaments;  (-  their  list  me  ■ 
from  diametrical  plane;  e - elevation  cf  free  vertices  above  tu:>lag*- 
axis;  h - distance  of  fuselage  axis  frcn  screen;  6 - radius  of 
cylindrical  fuselage,  work  [3]  snould  Le  used  as  a guide  in  selecting 
the  arrangement  of  the  free  vertex  filaments. 
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In  order  to  satisfy  the  condition  ct  ncnleakage  on  the  fuselage 
from  vortices  1 and  2 (Fig.  1)  in  the  case  cf  a wing  and  fuselage 
moving  in  an  infinite  fluid,  according  to  [/]  we  place  vortices  3 anu 
4 at  the  points  of  inversion  inside  the  fuselage.  In  order  to  satisfy 
the  boundary  condition  cf  ncnpenetr aticn  in  the  case  of  a system 
noving  near  a solid  screen,  we  map  on  it  ir  mirror  image  relative  tc 
the  plane  of  the  screen  the  preceding  vertex  system,  thus  obtaining  a 
hypothetical  vortex  system  (vortices  5-6).  In  order  to  eliminate  «-he 
viclated  condition  cf  nonpenetration  on  the  fuselage  while  preserving 
it  on  the  screen,  we  use  the  method  cf  successive  inversion  of 
vertices,  in  which  each  vertex  ci  the  hypothetical  system  (vcl*  x a0. 
Fig.  2)  is  mapped  in  relation  tc  the  surface  cf  the  fuselage,  placing 
at  the  point  of  inversion  vertex  Blc  which  is  mapped  relative  tc  the 
screen  - vortex  Ax.  Then,  we  perform  a second  approximation  - vertex 
A,  is  mapped  on  B2  and  Ea  - on  A?,  etc.  As  a result  of  this  chains  ct 
vortices  are  formed  inside  the  main  and  hypothetical  fuselages,  anu 
in  the  presence  of  vortex  A*0,  which  is  symmetrical  relative  the 
diametric  plane  of  vortex  A0,  the  same  chains  are  termed,  although 
this  time  symmetrical  in  relation  to  311  with  circulation  of  th- 
epposite  sign.  At  ter  the  y-erder  of  ap f to x i mat  i cn  the  "residual 
leakage"  or  the  fuselage  is  the  result  cf  vertices  % and  Aj , which, 
as  the  approximation  cider  increases,  move  closer  together,  reducing 
the  velocities  induced  by  and  thus  assurirg  nece  ss  air  y a c c u i a c y 

in  satisfying  th«-  ooundary  ccnai*ion  on  the  fuselage.  The  rcir;  hn 
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vortex  system  corresponding  to  the  studied  ^rctlem  is  shown  ir.  Fig. 
1r  where  dashed  lines  represent  the  chains  cf  vertices. 


lot  us  determine  vertical  downwash  ur  = at  point  M(x',  y', 

u© 

z')#  after  introducing  the  dimensionless  quantity 

r = R=y;  e=y;  h=j-,  x=j,  t,=j-  etc.  (1) 

After  applying  the  Bict-Savart  forirula  t.c  the  elements  or  the 
vertex  system  of  Fig.  1,  we  get: 


where 


V y"c*-cj+  (y-  i nf  i "Sr  ifi 


^ and  tn  are  the  coordinates  cf  vortices  1-d,  where  n is  the 
cf  the  vortex  according  tc  Fig.  1. 


The  values  of  the  coordinates  have  the  fern  of 


5,=  s.=  -Ct=-cs=?;  ?,=  ?*=e 

i + B 
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7*  and  are  the  coordinates  of  the  vortices  ir.  the  chain,  whe u x 
is  the  number  of  the  chain  (according  tc  Fig.  1)  ; y - the  number  or 
the  vertex  in  the  chain  (cr  order  of  ap pr oximat ion) . 

The  coordinates  of  the  vortices  ir.  the  chains  in  shifting  from 
approximation  y to  approximation  y + 1 vher  X = 6 and  X = 7 are 
calculated  according  to  the  forirulas: 


tr  = C*  . J*'_  _ _ HXiR+fj.) 

1 (<£)*+  (2E+?*)1  ’ > (g;y+  (2h+ ^7  ‘ 


(■>) 


Then,  the  coordinates  of  the  vortices  in  the  chains  for  the 
given  order  of  approximat  ion  (y  = ccnst)  are  related  to  ; c*  ; 
and  by  relationships: 


(h) 


The  coordinates  of  the  first  vortices  in  the  chains  (y  = 1)  when 
X = 6 and  x = 7 are  determined  ty  formulas  (5)  through  coordinates 
with  corresponding  indices  n = 6 and  n = 7,  which  are  calculated 
according  to  formula  (4). 
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Ey  using  the  Bernoulli  equation  and  allying  the  Biot-Savurt 
formula  to  the  elements  of  the  vortex  system  cf  Fig.  1,  we  get  a 
formula  tor  calculating  the  distribution  c£  tne  pressure  coefficient 
ever  the  surface  of  a cylindrical  fuselage: 


where 


fir  -yJ  i |/^a»e-^)+(Riine  - fy'i  {H) 

6 - is  the  polar  angle  of  the  studied  feint  on  the  fuselage  surface, 
reckoned  from  the  horizontal  plane. 


Ir.  calculating  quantities  and  f>  according  to  formulas  (2)  and 
(7)  we  should  limit  ourselves  to  approximation  •)  = k,  where  the 
increment  in  the  sum  with  respect  tc  y becomes  lower  than  the 
required  accuracy  of  calculating  quantities  ci  p. 


According  to  tnis  niethed  a calculation 
section  of  infinite  distance  downstream  for 
reduced  coefficient  of  pressure  over  the 


was  performed  for  a 
the  distribution  of  the 
surface  of  the  fuselage 
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(Fig.  3)  and  the  reduceu  vertical  uownwash  ~ along  the  winq  span 
(Fig.  4)  for  the  following  calculated  parameters: 

R = Y = 0,20;  ^ = ^-=0,25;  e = 0. 

Investigation  or  Velocity  Field  ori  FGDA  Device 

For  a section  of  infinite  distance  along  the  stream  the 
investigation  of  the  velocity  field  behind  a lew- flying  wing  with  a 
cylindrical  fuselage  can  oe" done  on  a plane  electrohydrody namic 
analogy  device  [the  EGD A ] , as  shewn  in  Fig.  5. 

The  symrcls  in  Fig.  5:  1 - sheet  cf  e lec t r cconduct ive  paper,  2 - 

tus,  3 - electrode  which  simulates  vortex,  4 - sheet  of 


elect  icconductive  paper  which  simulates  plane  undisturbed  tlow,  5 - 
direct  current  source,  6 - galvanometer,  7 - prebe. 

Simulation  of  the  free  vortex  flowing  eff  cf  the  wing  tic  is 
dene  ty  an  electrode  (source  or  electrical  field).  Observance  cf  t h- 
toundary  conditions  requires  that  the  ccntcur  (ct  fuselage  and 
screen)  in  this  analogy  be  made  from  a conductor.  The  symmetry  cf  the 
flow  in  relation  to  the  diametric  plane  is  assured  by  installing  a 
bus  in  this  plane.  In  this  case  the  fdee  lines  of  the  electrical 
field  will  correspond  tc  the  equ  ipoteri  t ia  1 lines  in  the  stream  cf 


D 


3 


8 


OC  - 05 


fluid,  the  eguipotential  lines  or  the  electrical  field  - tc  Kit  flew 
lines  in  the  fluid. 

The  posed  problem  for  a section  cf  inrinite  length  was  modeled 
cn  the  EGDA.  The  following  was  entained:  the  field  of  vertical 
downwashes  over  the  span  of  the  wing  (compared  with  calculations  cr. 
Fig.  4),  field  of  vertical  downwashes  at  location  of  tail  assembly 
(Fig.  b) , and  pattern  of  rlcw  lines  (Fic.  7)  fer  following 
parameters:  = 0,20;  K = 0,25  and  e * 0. 

CONCLUSIONS  1.  The  experiments  performed  cn  the 

elect rchydtodynamic  analogy  device  verified  the  proposed  method  ror 
calculating  the  velocity  field  behind  a low-flying  wing  wi+h 
cylindrical  fuselage  both  gual it  at i vely  and  quantitatively. 


2.  The  cylindrical  fuselage  substantially  effects  the  velocity 
field  near  the  fuselage,  within  a region  of  2-2.5  times  its  radius, 
decreasing  or  increasing  the  vertical  dewnwesh  depending  on  tne 
location  of  the  studied  [cint. 

3.  When  the  system  of  wing  and  fuselage  approach  the  screen,  the 


rarefaction  curve  on  the  surface  of  the  fuselage,  which  is  caused  h 
the  vertex  sheet  flowing  eff  of  the  wing,  decreases  on  one*  hand. 


EOC  =0531  9 

while  on  the  ot her  there  occurs  a redistr i t ut icn  in  this  curv-( 
result  of  wnich  additional  lifting  force  develops  on  the  fusel. 
This  deduction  must  also  qualitatively  le  verified  ty  the  fatt  . 
flew  lines  obtained  on  the  FODA. 
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